%
% studio dell'errore nel calcolo della derivata prima
% con le differenze in avanti e le differenze centrali
%

h = logspace(-16,0,1000);


% approssimazione della derivata del primo ordine

 df_1 = (exp(1+h)-exp(1))./h;

% approssimazione della derivata del secondo ordine

 df_2 = (exp(1+h)-exp(1-h))./(2*h);

% grafico dell'errore
 loglog(h,abs(exp(1)-df_1),'.b',h,abs(exp(1)-df_2),'.r');


